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摘要 :， 深 空 探测 中 探测 器 轨道 主要 采用 地 面 跟 踪 站 进行 测量 ， 无 线 电信 号 在 探测 器 与 地 面 跟踪 
站 间 传 输 时 会 受到 地 球 传播 介质 误差 的 影响 。 传 播 介 质 误差 主要 包括 地 球 的 中 性 大 气 时 延误 差 
和 电离 层 时 延误 差 等 。 针 对 传播 介质 误差 修正 在 深 空 探测 中 的 应 用 需求 ， 调 研 了 国内 外 地 基 测 量 
数据 传播 介质 误差 研究 进展 ， 重 点 阐述 了 深 空 探测 中 VLBI 地 基 测 量 数据 中 性 大 气 和 电离 层 时 
延 产生 的 机 理 和 修正 方法 。 针 对 关键 技术 及 未 来 发 展 趋势 进行 了 分 析 ， 以 提高 传播 介质 误差 修正 
精度 及 时 效 性 。 
x 88 dg: 深 空 探测 ， 中 性 大 气 时 延 ， 电 离 层 时 延 
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1 引 8 


至 20 世纪 50 年 代 末 ， 美 国 、 俄 罗斯 (前 苏联 )、 欧 盟 、 中 国 和 日 本 等 国家 和 组 织 已 经 
开展 了 250 多 次 深 空 探测 任务 。 探 测 目 标 基 本 涵盖 了 太阳 系 内 各 类 天 体 ， 有 些 探测 器 在 月 
球 、 火 星 和 小 行星 上 成 功 着 陆 或 返回 。 深 空 探测 工程 获得 了 大 量 科学 发 现 ， 为 人 类 认识 宇宙 
提供 了 重要 基础 ，。 

探测 器 发 射 升 空 后 ， 需 要 对 其 进行 跟踪 观测 。 美 国 采 用 国家 航空 航天 局 深 空 网 (deep 
space network, DSN) 跟踪 探测 器 ， 该 网 包括 3 座 经 度 间隔 120° 的 测控 站 和 1 个 位 于 加 利 
福 尼 亚 的 任务 操作 控制 中 心 组 成 。3 座 测控 站 是 美国 加 利 福 尼 亚 的 戈 德 斯 通 测控 站 、 西 班 牙 
境内 的 马德里 测控 站 以 及 澳大利亚 的 堪培拉 测控 站 。 这 种 覆盖 全 球 的 深 空 网 在 开展 探测 器 
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测控 时 随 着 地 球 的 转动 能 够 对 目标 保持 不 间断 的 跟踪 2 。 


长 基线 干涉 阵列 (very long baseline array, VLBA) 也 月 
测 土星 任务 时 采用 由 10 个 测 站 组 成 的 VLBA 阵列 进行 跟踪 “” 。 俄 罗斯 (前 苏联 ) 深 空 网 
个 测控 站 、2 个 指控 中 心 和 2 个 弹道 中 心 组 成 。3 个 地 面 站 分 别 是 乌 苏 里 斯 克 站 、 叶 夫 履 
托 里 亚 站 和 能 湖 站 ， 东 、 西 两 站 经 度 间隔 1009 左右 
天 操作 中 心 组 成 。3 个 测控 站 包括 澳 大 利 】 


3 


41 卷 


由 于 高 精度 的 测 角 精 度 ， 美 国 甚 


HiRes Wl as, Bd “KPa a” TR 


6] 


欧洲 深 空 网 由 3 个 测控 站 和 欧洲 航 


严 新 诺 舍 站 、 


西班牙 塞 布 里 洛斯 站 和 阿根廷 玛 拉 


库 站 "。 另 外 ， 欧 洲 其 长 基线 干涉 (very long baseline interferometry, VLBI) 网 (European 
日 本 宇宙 航空 研究 开发 机 构 管理 的 深 空 
网 主要 包括 白田、 谭 儿 岛 等 深 空 站 和 相模 原 控制 中 心 。 另 外 ,日 本 国立 天 文 台 的 4 个 VLBI 
测 站 也 参与 了 深 空 探测 ， 包 括 水 泽 、 入 来 、 小 笠原 和 石垣 岛 VLBI 站 ”“。 中 国 深 空 测控 系统 


VLBI network, EVN) 


主要 包括 3 个 测控 站 和 1 个 位 于 北京 的 控 外 


A 


VLBI 数据 处 理 中 心 位 于 上 海 ""。 
探测 器 飞行 过 程 中 采用 地 面 跟踪 站 进行 高 精度 跟踪 测量 以 准确 确定 其 飞行 轨迹 ， 尤 其 


IS 


南美 阿根廷 萨 帕 拉 "”"。 另 外 ， 中 国 


nM EX [4,5,8] 
也 参与 了 深 空 探测 任务 。 


Bras 3 个 测 


控 站 分 别 位 于 我 国 喀 什 、 佳 木 斯 


VLBI 网 (Chinese VLBI network, CVN) 同样 承担 了 
我 国 的 深 空 探测 任务 。4 个 射电 望远镜 跟踪 站 分 别 位 于 北京 、 上 海 、 昆 明和 乌鲁木齐 ，1 个 


在 中 途 轨 道 机 动 控 制 和 抵达 目标 星 后 


系 到 探测 任务 成 功 与 否 的 关键 。 
向 探测 器 的 距离 和 速度 ， 昨 
的 角度 信息 ， 从 而 完成 探测 器 的 定 轨 定 位 。 探 测 器 的 轨道 依赖 于 观测 量 的 动力 学 模型 及 测 站 
坐标 ， 影 响 测 距 、 测 速 的 误差 有 测 站 时 钟 的 不 稳定 性 、 传 播 介 质 误差 、 仪 器 误差 、 台 站 坐标 
自转 参数 引起 误差 等 。 长 弧 段 的 观测 带 来 测量 误差 的 积累 ， 影 响 测量 精度 。 随 着 


误差 和 地 球 


深 空 探测 任务 的 不 断 深入 ， 尤 其 在 探 
信息 。 


20 世纪 80 年 代 ， 来 自 射 


制 动 等 关键 弧 段 。 探 测 器 能 否 
早期 的 深 空 探测 器 定 轨 采 月 


HIE 


Ne 


顺利 到 达 指 定位 置 是 关 
测速 技术 测量 位 于 视线 方 


[UT 
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通过 长 时 间 长 弧 段 的 多 普 勒 速度 测量 ， 间 接 获 取 探测 器 在 平面 上 


测 器 着 各 


和 返回 等 关键 任务 段 需要 更 精确 的 探测 器 位 置 


HA SOMA VLBI 技术 开始 应 用 于 深 空 探测 领域 。 该 技术 由 


两 个 或 两 个 以 上 的 天 线 分 别 在 同一 时 刻 接收 同一 探测 器 或 河 外 射电 源 的 信号 ， 然 后 进行 相 


关 处 理 求 出 观测 值 。 
PER, 


REME HBR ELS RA AN HE 


频 系统 不 稳定 性 、 
观测 探测 器 附近 的 河 外 射 : 


这 种 干涉 测量 的 优点 是 基线 长 度 可 达 几 和 王 干 米 ， 极 大 地 提高 了 空间 角 分 
可 以 更 精确 测定 探测 器 角度 信息 。VLBI 技术 的 主要 测量 误差 分 为 两 类 : 一 类 是 由 于 


角 引 起 的 误差 ， 称 为 模型 与 参数 误差 ， 男 一 类 是 测 站 时 
测 站 信号 接收 系统 引起 的 信号 相位 抖动 和 传播 介质 时 延 影响 等 。 通 过 交 蔡 
电源 ， 可 以 消除 共同 的 仪器 和 测 站 误差 等 。 由 于 传播 介质 引起 的 时 


延误 差 随时 间 而 变化 ， 并 且 探 测 器 和 射电 源 间 的 角 距 有 数 度 至 十 几 度 ， 导 致 不 同 路 径 上 介质 
时 延 差 异 ， 见 图 1。 


影响 深 空 探测 信号 的 传播 介质 包括 地 球 空间 环境 和 
深 空 探测 地 基 测 量 数据 中 地 球 空间 环 


f 


IN 


延 修正 方法 。 


星际 空间 环境 。 本 文 主 要 研究 了 


iB] rp M 


fT 
E 大 气 和 电离 


:时 延 产 生 的 机 理 、 修 正方 法 及 


究 现 状 ， 同 时 分 析 了 时 延 修正 关键 技术 及 未 来 发 展 趋势 。 


As 


简要 介绍 了 太阳 等 离子 体 引 起 的 
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射电 源 e 


1 无线 电信 号 传播 示意 图 


2 HEKA 


本 文 所 指 中 性 大 气 主要 是 距离 地 球 表面 约 50 km 的 大 气 ， 包 括 对 流 层 和 平流 层 。 对 流 
层 的 高 度 因 纬度 而 不 同 ， 低 纬度 地 区 为 17 ~ 18 km， 中 纬度 地 区 为 10 ~ 12 km， 极 地 为 8 
~ 9 kmo 
2.1 ”中 性 大 气 时 延误 差 产 生机 理 

大 气 以 两 种 方式 影响 电磁 波 的 传播 : (1) 与 在 真空 中 相 比 电磁 波 的 传播 速度 减 慢 ，(2) 
传播 路 径 由 直线 变 成 弯曲 。 这 两 种 影响 都 是 传播 路 径 上 大 气 折 射 指数 变化 造成 的 。 信 和 号 传播 
在 时 间 上 的 延迟 等 效 于 传播 路 径 长 度 的 增加 ， 增 加 的 路 径 表示 为 : 


AL-aAt-G-« | -a= /nods-G d (1) 
LC L 


HP, AL 为 由 于 中 性 大 气 时 延 增加 的 路 径 ，co 为 真空 中 的 光速 ，At 为 信号 实际 传播 时 间 ， 
G 是 探测 器 和 接收 端的 几何 路 线 长 度 ， 工 为 弯曲 路 径 ，e 为 信号 在 大 气 中 的 传播 速度 ，n(s) 
是 大 气 折射 指数 。 


AL= [ ms) 1145-6 /pa-lde+s-G ; (2) 


S 为 沿 工 的 路 径 长 度 。 右 侧 第 一 项 积分 项 为 减 慢 项 ， 第 二 项 是 弯曲 项 ， 为 弯曲 路 径 与 直线 
路 径 之 差 。 弯 曲 延迟 在 探测 器 高 度 角 较 小 (< 5°) 时 会 有 较 大 差 值 ( 约 2 m)， 但 相对 于 传播 
速度 延迟 来 说 要 小 得 多 。 当 探测 器 高 度 角 为 15° 时 ， 有 十 几 厘米 水 平 ， 并 随 着 高 度 角 的 增 
大 ， 弯 曲 延迟 趋 近 于 0， 可 忽略 不 计 "。 
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H 
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性 大 气 时 延 的 大 小 主要 取决 于 电磁 波 传播 路 径 上 大 气 的 厚度 以 及 各 点 的 气象 参数 ( 包 
括 温 度 、 湿 度 和 压强 )。 具 有 以 下 特点 : (1) 时 延 变化 幅度 较 大 ， 在 天 项 方向 约 2.3 m， 低 高 
度 角 时 可 达到 30 m; (2) 在 30 GHz 以 下 具有 非 色散 性 所， 不 能 利用 频率 组 合 的 方式 消除 ; 
(3) 大 气 中 水 汽 分 布 随机 性 强 ， 难 以 建 模 ，(4) 不 同 地 区 、 不 同 季节 差异 明显 ， 全 球 建 模 难 
度 大 。 
中 性 大 气 时 延 随 着 高 度 角 而 变化 ， 可 近似 为 天 项 时 延 与 观测 高 度 角 正弦 倒数 的 乘积 。 基 
T VLBI 技术 对 探测 器 测控 时 采用 差分 观测 ， 探 测 器 与 河 外 射电 源 之 间 存 在 一 定 的 角 距 。 所 
以 ， 中 性 大 气 时 延 与 高 度 角 、 探 测 器 河 外 射电 源 角 距 高 度 相 关 ”。 如 下 式 所 示 : 


E 1 1 

ATA clsinelp  sinelc| ' (3) 
AP, Amas 是 不 同 观 测 高 度 引 起 的 较 差 时 延 ，s 为 大 气 天 顶 时 延 的 精度 ，elr 和 elo 分 别 为 
探测 器 和 河 外 射电 源 的 观测 高 度 角 。 
图 2 给 出 了 中 性 大 气 时 延误 差 对 定 轨 时 延 的 影响 随 高 度 角 和 探测 器 与 河 外 射电 源 角 距 
的 变化 关系 。 图 上 中 性 大 气 天 顶 时 延误 差 设 为 2 ecm， 和 角 距 分 别 设置 为 22、5° 和 8°*。 当 高 
度 角 为 10"， 角 距 为 8" 时 ，2 cm 的 中 性 大 气 天 顶 时 延误 差 对 定 轨 时 延 影 响 会 达到 170 ps 
左右 。 


时 延误 差 /ps 
I 


80r 
40r 
0 1 L 1 Sas 
10 20 30 40 50 60 70 
高 度 角 /(?) 


ik: 红色 代表 角 距 8”; 黑色 代表 角 距 5" ， 蓝 色 代表 角 距 2。 


2 ”中 性 大 气 时 延误 差 随 高 度 角 和 角 距 的 影响 


所 以 ， 选 择 小 角 距 的 参考 源 ， 避 免 在 低 高 度 角 观测 ， 可 以 降低 中 性 大 气 时 延 影 响 。 然 
而 ， 小 角 距 条 件 下 往往 无 法 找到 满足 强度 要 求 的 射电 源 。 为 了 增 大 基线 长 度 ， 测 站 间 往 往 相 
阳 几 千 千 米 ， 在 同时 观测 的 情况 下 不 能 达到 所 有 测 站 与 探测 器 间 高 度 角 符合 要 求 。 另 外 ， 中 
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性 大 气 时 延误 差 与 测 站 坐标 、 接 收 机 钟 差 高 度 相 关 ， 因 此 不 能 采用 参数 估计 方法 消除 中 性 大 
气 时 延误 差 的 影响 “ ， 使 用 外 部 数据 进行 中 性 大 气 时 延 修 正 是 必须 的 校准 过 程 。 
2.2 ”中 性 大 气 时 延 修 正法 

本 文 对 国内 外 测控 站 中 性 大 气 时 延 修 正方 法 进行 调研 分 析 ， 目 前 常用 的 方法 包括 天 顶 
时 延 十 映射 函数 法 和 基于 微波 辐射 计 (Water Vapor Radiometer, WVR) 修正 法 等 。 
2.2.1 中 性 大 气 天 顶 时 延 十 映射 函数 法 

将 探测 器 方向 上 的 信号 时 延 表 示 为 天 顶 时 延 和 与 高 度 角 有 关 的 映射 函数 的 组 合 ， 是 目 
前 深 空 探测 中 地 基 测 量 数据 中 性 大 气 时 延 有 效 的 修正 方法 。 
2.2.1.1 中 性 大 气 天 顶 时 延 

(1) 美国 探测 器 中 性 大 气 天 顶 时 延 修正 法 

20 世纪 70 年 代 开 始 ， 美 国 深 空 网 中 性 大 气 天 顶 时 延 修正 方法 主要 采用 Chao 提出 的 季 
节 模型 修正 “…。 与 探 空气 球 结果 比较 ， 该 模型 天 顶 湿 时 延 1o 精度 约 2 cm。 由 于 不 需要 实 
测 资料 ， 该 模型 可 作为 预报 模型 应 用 于 实时 中 性 大 气 时 延 修正 ， 

随 着 探测 器 定 轨 精 度 的 不 断 提 升 ， 原 有 的 季节 模型 已 无 法 满足 中 性 大 气 时 延 修 正 的 要 
求 。 美 国 喷气 推进 实验 室 (Jet Propulsion Laboratory, JPL) 提出 中 性 大 气 时 延 修正 多 项 式 
模型 "”"， 采 用 时 间 归 一 化 傅 里 叶 序 列表 示 干 时 延 (zenith hydrostatic delay, ZHD) 和 湿 时 延 
(zenith wet delay, ZWD)， 如 式 (4) Aras: 


Fs(X)= Ao + Ai: cos(X) + By: sin(X) + A2: cos(2X) + B5: sin2X) +--+ 
An -cos(NX)+ By -sin(NX) , 


Xon (=£) | (4) 


Jh, A, B; WAHAB N 是 傅 里 叶 序 列 展开 式 的 阶 数 ，7 RINT, S 是 开 
始 时 间 ，P 是 基本 模式 周期 ， 单 位 都 为 s。 

其 中 ， 多 项 式 模型 的 一 阶 项 是 Chao 提出 的 季节 模型 “ ”， 参 数 以 时 间 归 一 化 傅 里 叶 序 
列表 示 。 为 避免 双重 校准 ， 在 季节 模型 校准 的 基础 上 ， 采 用 近 实 时 气象 资料 和 全 球 定位 系统 
(Global Positioning System, GPS) 观测 数据 计算 中 性 大 气 天 顶 时 延 、 干 时 延 和 湿 时 延 的 增 
量 校正 时 间 序 列 ， 并 采用 时 间 归 一 化 多 项 式 进行 表示 ， 如 式 (5) 所 示 : 


Py(X) = Cot C; X c9 X cc Cy- XP, 


x- (2-5). (5) 


HH, C 是 拟 合 多 项 式 系数 ，N 是 多 项 式 展开 式 的 阶 数 ，T 是 标 校 时 间 ，S 是 开始 时 间 ， 
是 结束 时 间 。 
全 球 导航 卫星 系统 (Global Navigation Satellite System, GNSS) 和 中 性 大 气 折 射 时 延 紧 
密 相 关 ， 一 方面 卫星 信号 在 穿 过 中 性 大 气 层 时 会 受到 其 折射 的 影响 带 来 一 定 的 时 延误 差 ， 另 
一 方面 ， 利 用 该 影响 可 通过 GNSS 技术 反 演 获得 中 性 大 气 时 延 量 。 由 于 GNSS 技术 具有 全 


ies 
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天 候 、 高 采样 、 高 精度 和 不 受 天 气 条 件 影响 等 优点 ， 因 此 它 是 中 性 大 气 时 延 改正 的 重要 手 


段 。 


该 技术 成 功 应 用 到 深 空 探测 地 基 测 量 数 据 中 愧 


(slant path delay, SPD) 可 以 通过 天 顶 时 延 、 映 射 函数 和 梯度 进行 表示 ， 


SPD = mfr: ZHD +mfu: ZWD 4t mf,-(Gycoso-- Ggsino) ， 


大 气 时 延 修正 。 传 播 路 径 上 中 性 大 气 时 延 


LA PAX: 
(6) 


HB. mf, 和 mf 分 别 为 干 时 延 和 湿 时 延 映射 函数 ，2ZHD 和 ZW D 分 别 为 天 顶 方 向 干 时 


AERE 
位 角 。 


美国 深 空 网 基于 GPS 技术 反 演 中 性 大 气 时 延 的 方法 是 ; 


参数 ， 采 用 最 小 二 乘法 或 卡尔 曼 滤 波 法 ; 
的 GIPSY-OASIS。 该 软件 与 伯 
平均 差异 为 6 mm， 标 准 差 为 3.2 mm. | 


为 了 得 到 湿 时 延 ， 基 于 3 


尔 尼 大 学 


行事 后 或 实时 处 到 


TE, mf, 是 梯度 映射 函数 ，Gx 和 Ge 分 别 是 南北 向 和 东西 


各 中 性 大 


ES 


T 


向 水 平 梯度 ，a 是 方 


气 天 项 总 时 延 作为 待 估 


B. 数据 处 理 软件 采用 JPL 研发 


研发 的 Bernese 软件 估计 的 
于 采用 GNSS 技术 只 能 获取 中 性 大 气 天 项 总 时 延 
置 站 气象 数据 采用 Saastamoninen 模型 计算 干 时 延 ，GNSS 观测 数 


据 反 演 的 中 性 大 气 总 时 延 扣除 干 时 延 后 得 到 深 空 站 湿 时 延 信 息 。 


2001 年 美国 火星 探测 器 “奥德赛 ” 定 轨 时 ， 地 基 测 量 数据 中 性 大 气 时 延 由 
中 性 大 气 时 延 采 用 上 段 介 
获取 的 中 性 大 气 时 延 对 深 空 站 ADOR、VLBI 和 VLBA 观 测 值 进 行 校准 。 
12h ZTD 方法 ， 可 应 用 于 深 空 


出 了 可 用 于 延迟 1 和 


ZTD 产品 修正 精度 在 厘米 量 级 ， 


esl 


用 到 “ 
KI 


FE 西 尼 ”探测 器 轨道 确定 


， 美 国 VLBA 也 用 了 


的 GPS 技术 处 到 


延迟 12 h 提供 的 ZTD 产品 修了 
20] 


探测 器 定 轨 。 美 国 
个 强 射 电源 标 校 中 性 大 气 天 顶 时 延 。 该 方法 通过 在 大 范围 


方法 修正 。2007 4 


VLBA 测定 “ 


W 


中 性 大 气 天 项 湿 时 延 的 


[18] 


JPL 提供 ， 


, JPL 采用 上 述 GPS 技术 


同年 ，JPL 提 


AAI BE IE, HEIR 1 h 提供 的 
EEF lcm. 该 方法 应 


FE 西 尼 ” 探 测 器 时 采用 约 15 


高 度 角 区 间 快 速 观 测 强 射电 源 以 


获取 多 波段 延迟 ， 然 后 通过 Chao 映射 函数 投影 到 探测 器 观测 方向 “2 。 


(2) 欧洲 探测 器 中 性 


欧洲 深 空 网 常用 的 中 改 


气象 资料 计算 中 性 大 气 时 延 ， 可 用 于 事后 修正 。 由 于 大 气 中 水 蒸气 的 


EX 


E FUR TIS] ET TEV 
秆 延 修正 方法 是 Saastamoinen WW, 


该 方法 基于 实测 
非 均 匀 分 布 ， 该 方法 获 


BU" ZHD 精度 较 高 ， 但 ZWD 精度 较 差 ， 在 厘米 量 级 。 基 于 合作 协议 ， 欧 洲 深 空 网 中 性 大 
气 时 延 修正 也 由 美国 JPL 提供 支持 。 


近年 来 ， 为 了 不 依赖 外 部 数据 ， 基 
一 个 独立 的 中 性 大 气 时 延 修 正平 台 。 


于 JPL 提出 的 多 项 式 修正 算法 ， 欧 洲 深 空 网 开发 了 
主要 内 容 包 括 : (1) 导航 办 公 室 基于 


GNSS 技术 反 演 提 


供 深 空 站 48 h 时 段 长 度 的 ZTD 和 ZWD. (2) 每 个 欧洲 深 空 站 配置 气象 仪 ， 提 供 1 min 采 


样 的 温度 、 


湿度 和 压强 信息 ， 基 于 Saastamoinen 模型 计 
全 球 压 强 温 度 模 型 (global pressure and temperature, GPT) 计 


ZHD。 气 


压 数 据 不 可 用 时 ， 采 用 


3 (3) 改正 由 于 深 空 站 


天 线 和 国际 GNSS 服务 组 织 (International GNSS Service, IGS) 跟踪 站 高 度 差 引 起 的 ZHD 


差异 ， 深 空 站 的 参考 点 为 天 线 方位 轴 


与 仰角 轴 相 交 的 点 。 其 体 方法 参考 文献 [17]。(4) 针对 


参与 解 算 的 GNSS 测 站 数 不 同 ， 分 快速 模式 和 最 终 模 式 。 最 终 模 式 延迟 4d 计算 。(5) 中 性 


大 气 时 延 多 项 式 阶 数 为 4 阶 ， 采 用 最 小 二 乘 估计 多 项 式 系数 。 图 3 是 


欧洲 导航 办 公 室 开发 的 
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中 性 大 气 时 延 修正 平台 计算 ZWD, ZHD 和 标 校 中 性 大 气 时 延 的 时 间 方 案 ， 每 天 14:30(UT) 
启动 计算 当天 12:00(UT) 往 前 48 h 观测 数据 。 


HX 14:30 
E  : 


[13] 


3 ”欧洲 导航 办 公 室 计算 ZWD、ZHD 和 标 校 中 性 大 气 时 延 的 时 间 方 案 


欧洲 EVN 也 用 于 探测 器 轨道 跟踪 ， 比 如 “织女 哈雷 彗星 ”探测 器 探测 金星 大 气 时 
期 ,“ 惠 更 斯 ”探测 器 降落 在 土 卫 六 表面 期 间 ””， 欧 空 局 “Smart-1” 探 测 器 在 月 球 表面 控 
制 期 间 ， 美 国 火星 探测 车 B 最 后 巡航 阶段 ”,“ 卡 西 尼 ” 号 土星 探测 器 “， 以 及 最 近 的 “人 金 
ERE” RW” KA EVN 进行 轨道 确定 时 ， 中 性 大 气 时 延误 差 有 两 种 修正 方法 : (1) 
采用 经 验 模 型 和 维也纳 映射 函数 (Vienna mapping functions 1, VMF1); (2) 采用 射线 追踪 
技术 ”" ， 该 技术 在 后 面 进行 具体 介绍 。 在 高 度 角 较 低 或 探测 器 射电 源 角 距 较 大 时 采用 第 
种 方法 ， 例 如 “金星 快车 ”跟踪 任务 EM081lcD 。 

(3) 日 本 探测 器 中 性 大 气 天 顶 时 延 修 正法 

2007 £, HÆ SELENE 探测 器 发 射 升 空 ， 日 本 射电 天 体 测量 VLBI 网 和 国际 VLBI 
网 (包括 我 国 上 海 、 乌 鲁 木 齐 ， 澳 大 利 亚 霍 巴特 和 德国 韦 茨 拉 尔 ) 采用 多 频 VLBI 技术 进 
行 SELENE 探测 器 轨道 确定 。Liu 等 人 基于 VLBI 同 波 速 技术 进行 SELENE 探测 器 的 轨 
道 确定 时 ， 中 性 大 气 时 延 采 用 GPS 技术 修正 。 其 中 ，ZHD 采用 测 站 气压 值 计算 ， 精 度 为 
1mm; ZTD 采用 GpsTools 软件 处 理 ， 数 据 处 理 时 的 轨道 和 钟 差 文件 采用 IGS 网 站 最 终 
产品 ，ZTD 精度 为 10 mm™ ”。 另 外 ， 基 于 合作 协议 ， 日 本 探测 器 也 由 美国 深 空 站 跟踪 
测 轨 。 

(4) 我 国 探测 器 中 性 大 气 天 顶 时 延 修 正法 

我 国 深 空 探测 与 美国 和 欧洲 相 比 起 步 较 晚 ， 但 发 展 速 度 快 。 针 对 中 性 大 气 时 延 修正 ， 
2014 年 ， 赵 军 等 人 基于 测 站 历史 气象 数据 拟 合 折射 衰减 系数 ， 建 立 我 国 测控 站 上 空中 性 大 
气 折射 率 模型 ， 并 利用 实测 地 面 折 射 率 对 统计 模型 进行 优化 ， 该 模型 成 功 应 用 于 探测 器 长 管 
外 测 数 据 的 实时 修正 ”。2019 年 ， 北 京 航天 飞行 控制 中 心 构建 了 一 种 适用 于 深 空 探测 实时 
任务 下 的 区 域 高 精度 中 性 大 气 时 延 模型 。 首 先 ， 通 过 中 性 大 气 时 延 实测 值 与 模型 值 的 比 对 
确定 Saastamoinen 模型 天 顶 时 延 估 计 的 修正 系数 ， 修 正 后 中 性 大 气 天 顶 时 延 偏差 由 1.24 ns 
降 至 约 0.9 ps; 其 次 ， 为 突破 Saastamoinen 模型 气象 数据 约束 ， 基 于 高 程 修正 的 UNB3m 
模型 ， 建 立 了 测 站 区 域 大 气 参数 估计 模型 ， 大 气压 强 、 温 度 、 水 汽 压 的 估计 偏差 分 别 约 为 
390 Pa, 6.7 K, 63 Pa， 由 大 气 参 数 模型 值 与 实测 值得 到 的 中 性 大 气 天 顶 时 延 偏 差 约 为 14 ps。 


| 
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采用 该 方法 ,“ 娣 娥 四 号 ”探测 器 实时 定 轨 残 差 与 事后 中 性 大 气 时 延 实测 值 条 件 下 的 结果 基 
A 
RE CVN 也 对 深 空 探 测 器 进行 跟踪 观测 ， 比 如 “嫦娥 ”系列 探测 器 和 “天 问 一 号 ” 火 
星 探测 器 。 近 年 来 ， 基 于 北京 、 上 海 、 昆 明和 乌鲁木齐 测 站 3 年 ZTD 序列 ， 中 国 科学 院 上 
海天 文 台 研 发 了 CVN 测 站 本 地 化 天 项 时 延 模型 “” ， 该 模型 不 需要 实测 数据 ， 可 以 应 用 于 
VLBI 观测 传播 介质 误差 实时 修正 。 基 于 “嫦娥 五 号 ”实测 数据 分 析 ， 该 模型 计算 的 中 性 大 
气 天 顶 总 时 延 在 4 个 测 站 上 的 均 方 根 (root mean squire, RMS) 分 别 为 0.11, 0.14, 0.03 和 
0.08 ns。 采 用 该 模型 后 低 高 度 角 (10° ~ 28°) 时 VLBI 时 延 残 差 约 0.62 ns。 模 型 表达 式 为 : 


j 


tho . tho 
ZTD = ZTD + A cos ( hoy 2) + Bisin ( noy 2r) 4 


365.25 x 24 365.25 x 24 
tho . tho 
A Y_ag) +B Y_ag ) 4 
oo, E x 24 ) Xn m x 24 ) 
tho š tho 
A y B se EY eet 
UNE (ss x sao") TU (a x z0r) i (7) 


Kb, ZTD 为 测 站 上 中 性 大 气 天 顶 时 延 ，2ZTDo， A Bi 为 单个 站 点 上 的 拟 合 系数 ， 刀 ay 是 
年 内 小 时 ， 使 用 时 需要 根据 测 站 确定 对 应 参数 及 对 应 时 间 点 。 

为 了 进一步 提高 中 性 大 气 时 延 修 正 精度 ， 中 国 科学 院 上 海天 文 台 基于 GNSS 实测 数据 
实现 了 1 更 新 的 GNSS 反 演 中 性 大 气 时 延 修正 技术 ， 并 应 用 于 VLBI 技术 定 轨 事 后 3 h 模 
式 。 通 过 “ 娣 娥 五 号 ”实测 数据 分 析 ， 在 低 高 度 角 (10° ~ 28°) 观 测 时 VLBI 时 延 残 差 RMS 
可 达到 0.26 ns 
2.2.1.2 ”映射 函数 

20 世纪 70 年 代 ， 美 国 深 空 网 采用 Chao 等 人 提出 的 映射 函数 模型 ， 将 中 性 大 气 天 顶 时 
延 投影 到 探测 器 观测 方向 77, 20 世纪 80 年 代 Lanyi, Davis 和 Herring 分 别提 出 了 半 解 
析 映 射 函 数 、CFA2.25 8 MTT "Begg. 20 世纪 90 年 代 ，Neill 提出 了 NMF ”， 严 豪 
健 等 人 "新 建 了 UNSW 模型 。 对 于 CFA.2.2 和 UNSW 模型 ， 当 地 面 气象 记录 偏离 模型 规 
范 值 较 大 时 ， 映 射 函 数 精度 会 很 快 下 降 。MTT 和 NMEF 模型 可 以 较 好 地 描述 大 气 变化 范围 
的 大 气 剖 面 的 映射 函数 。 但 由 于 探 空气 球 的 局 限 性 ， 该 模型 在 全 球 不 同 地 区 修正 精度 差异 
较 大 。 

美国 深 空 网 、 欧 洲 深 空 网 、 日 本 VLBI WARE CVN 主要 采用 NMF 模型 将 ZTD S 
影 到 探测 器 观测 方向 “”。 欧 洲 EVN 采用 VME1 模型 获取 探测 器 观测 方向 中 性 大 气 时 延 
修正 量 "。 北 京 航天 飞行 控制 中 心 通过 分 析 NMP 映射 函数 规律 ， 确 定 干 映射 参数 a 为 响应 
的 主要 影响 参数 ， 进 一 步 通 过 最 小 二 乘法 计算 出 参数 a 的 最 佳 值 。 结 果 表 明 在 高 度 角 为 10° 
时 ， 中 性 大 气 时 延 偏差 约 0.3 ns， 相 对 缩小 近 1 个 量 级 。 该 方法 应 用 于 “嫦娥 四 号 ”实时 任 
务 ”。 近 两 年 ， 我 国 CVN 采用 GMP 模型 进行 投影 ”。 

21 世纪 以 来 ， 随 着 数值 天 气 模型 的 建立 ， 维 也 纳 工业 大 学 发 展 了 VME 系列 模型 “”。 
该 系列 模型 包括 第 一 代 (GMF/VMF1/VMF1-FC)， 第 二 代 (GPT2/GPT2w/VMF2) 和 第 三 
ÍR (GPT3/VMF3/VMF3-FC)， 第 三 代 模 型 精度 最 好 。 其 中 VMF1, VMF2 和 VMF3 需要 


Lm 
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实测 格 网 数据 ， 可 
实时 修正 ，VMF1-FC 和 
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以 在 事后 应 用 ;GMF, GPT2, GPT2w 和 GPT3 是 经 验 模型 ， 
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可 以 用 于 


VMF3-FC 为 采用 格 网 预报 产品 构建 的 模型 ， 同 样 可 以 应 用 于 实时 


修正 ， 但 该 类 格 网 预报 产品 需要 权限 
实时 中 性 大 气 映射 函数 GPT3、VMF3-FC 和 站 处 VMF3-FC 的 建 模 精度 进行 了 评估 。 
VMF3-FC 精度 最 优 ; VMF3-FC 次 之 ， 在 无 站 处 产品 的 情况 下 推荐 使 用 ; 
最 差 ， 但 其 可 内 内 在 软件 中 ， 在 不 能 下 载 格 网 预报 产品 的 情况 下 寺 
考虑 采用 该 类 模型 提高 探测 器 观测 方向 中 必 
基于 WVR 的 中 性 大 气 时 延 修正 方法 

目前 最 精确 的 中 性 大 气 时 延 修 正方 法 是 基于 


空 站 可 


2.2.2 


才 可 以 获取 。 武 汉 大 学 采用 射线 追踪 法 对 最 新 一 代 模 型 


站 处 
GPT3 整体 表现 
PEN. ÆT, R 


观测 方向 上 的 亮 温 ， 然 后 通过 专用 反 演 算法 


的 压强 基于 反 演 算法 获取 。 


EXC] GE 


BEE. 


WVR 实测 ， 该 方法 可 以 直接 测量 探测 器 
估计 中 性 大 气 湿 时 延 。 干 时 延 通过 WVR 观测 


失效 。 这 也 是 该 方法 应 用 在 深 空 探测 的 一 个 严重 限制 因 


美国 JPL 最 早 发 展 了 WVR。 


该 技术 受 液态 水 的 影响 ， 在 大 雨 或 多 云天 气 条 件 下 反 演 算法 可 能 
素 。 


1995 年 ，Keihm 对 加 利 福 尼 亚 戈 德 斯 通 站 一 年 的 WVR 


数据 与 本 站 预报 模型 进行 统计 分 析 ， 中 性 大 气 修正 精度 在 1 h 内 的 时 间 尺 度 上 可 以 提高 


30%， 超 过 1 h 则 与 预报 模型 精度 相当 ”。 


WVR， 安 装 于 深 空 站 ， 用 于 “ 卡 西 


1000 ~ 10000 s 时 间 尺 度 上 的 亮 


0.05 K, 1a 时 间 尺 度 上 的 亮 温 稳 定性 为 0.15 K. 
DSN 各 站 上 配备 了 IPL E 
度 校 准 中 性 大 气 波 动 。 


完成 深 空 探测 任务 ， 美 国 
100 ~ 10000 s 的 时 间 尺 度 上 高 精 


2003 年 ，JPL 的 Tanner 和 Riley 设计 了 新 型 
尼 ” 探 测 器 定 轨 中 校准 中 性 大 气 时 延 的 影响 。 通 过 两 
台 辐 射 计 的 相互 比较 和 近 1 a 的 连续 曲线 标定 ， 对 辐射 计 的 性 能 进行 了 评估 。 该 WVR 在 
温 稳定 性 优 于 0.01 K, 30 d 时 间 尺 度 上 的 亮 温 稳定 性 为 
2007 E, NTL “FAE” RI 


BCH) 7A RR an fa EK WVR, HF 


实验 表明 ， 在 短 基 线 干 涉 测量 时 间 


尺度 2000 ~ 10000 s 中 ， 湿 时 延 校准 的 相位 残 差 阿 伦 标 准 差 优 于 2 x 10719. ZEAE IRR E 
下 ，WVR 可 以 将 水 蒸 气 引起 的 相对 时 延 校准 到 1 mm (RMS) 以 上 ， 绝 对 时 延 校准 到 3 mm 


以 上 SW. 
欧洲 深 空 站 各 配置 了 2 套 德国 


2014 年 11 月 安装 在 
干扰 最 小 。 该 WVR 由 2 个 


6 班 牙 站 ，WVR E 


水 汽 吸 收 波段 ， 另 7 个 通道 


研发 的 WVR (型 号 RPG-HATPRO)。 
E 离 深 空 站 射电 望远镜 天 线 100 m， 确 保 与 主 天 线 的 
接收 机 共 14 个 通道 组 成 ， 其 中 7 个 通道 
于 测量 50 ~ 60 GHz 氧气 吸收 波段 。 这 
供 高 热 稳定 性 和 测量 精度 。 通 过 对 安装 在 荷兰 卡 鲍 的 WVR (型 号 RPG-ATPROP) i 


第 1 套 WVR 于 


于 测量 20 ~ 30 GHz 
Yi n] USE 2f hp 


究 ， 发 现 天 线 波 速 对 中 性 大 气 时 延 修正 精度 的 影响 较 小 ， 射 电 天 线 与 WVR 最 佳 相对 距离 
为 25 m， 两 个 天 线 的 指向 偏差 为 毫米 量 级 ”。 


[46] 


2014 年 Graziani 等 人 


减少 到 其 零点 偏 移 值 的 5596. cd 


在 探测 器 定 轨 中 校准 中 性 大 气 变化 的 作用 。 
站 天 线 波束 相交 ， 风 量 中 各 点 之 间 的 虽 


(相距 10 km) 均 具 有 较 好 的 一 致 特 
提供 技术 参考 。 


E", JR ERO 


采用 阿 伦 标准 差 (Allan standard deviation, ASD) 评估 了 WVR 
通过 为 WVR 波束 引入 指 
E 离 减 小 ， 从 而 减少 残余 ASD。 测 试 表明 产生 的 ASD 
RZY WVR 与 ITU-R 模型 "和 附近 美 
本 联合 打造 的 “ 贝 皮 科 伦 布 号 ”水 星 探 测 


可 偏 移 ， 使 其 与 深 空 


玉 | 


深 空 站 WVR 
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2013 年 左右 ， 中 国 
WVR. 2020 年 ， 上 海 VLBI 站 本 


LA- EE 


喀什 和 佳木斯 深 空 站 分 别 安装 了 


天 文学 进展 


E 


41 卷 


7 


d 


H 


电波 传播 研究 所 生产 的 


jt VLBI 站 分 别 配 置 一 套 中 国 
精度 可 以 达到 0.2 K^ 


Be Fe RED FE TA 


10 mm， 中 国电 波 传播 研究 所 HZD-X Ili] HEE 1 
2.2.3 其 他 方法 
探 空 资料 


。 通 过 与 GNSS 技术 对 比 ， 德 国 RPG-G5 天 顶 湿 时 延 人 
F 精 度 优 于 15 mm. 


| 和 数值 天 气 模型 (numerical weather model, NWM) 资料 等 可 以 提供 


WVR (型 号 RPG-G5)， 北 京 、 昆 明和 乌 鲁 
E 产 的 WVR (型 


E 


2-5 HZD-X)， 亮 温 长 期 测 
Ef BE 


P3 


三 维 气象 


信息 ， 为 计算 大 气 的 折射 率 场 提供 了 可 能 ， 利 用 射线 追踪 技术 在 折射 率 场 的 基础 上 对 信和 号 


传播 路 径 上 的 折射 率 进行 积分 ， 可 
可 以 提供 数值 天 气 预 报 资料 ， 如 
预报 中 心 发 布 的 气候 预报 系统 上 


保证 。 
约 16 mm 


p [8,27] 


方法 


ESI 
[26] 


， 在 EVN 探测 “金星 ”任务 


2016 年 ， 意 


> sy 


大 利 学 者 提 


以 直接 获得 大 气 对 
欧洲 中 期 天 气 预 报 中 ， 
了 分 析 资 料 等 。 当 前 的 NWM 模型 中 将 更 多 小 尺度 天 气 现象 
融入 模型 的 运行 ， 极 端 天 气 情况 下 通过 射线 追踪 技术 计 
射线 追踪 技术 获取 的 ZTD 与 GPS 技术 反 演 ZTD 对 比 ， 偏 差 为 0.9 mm, RMS 
FP， 高 度 角 较 低 或 探测 器 射电 源 角 距 较 大 


出 天 气 预报 数值 模型 与 微 物理 


中 


Hc HIT RE 
Vy UR 


1 
TES 


FE 
Ea 


目前 多 个 
下 分 析 资 料 ， 美 


家 和 地 区 
到 家 环境 


大 气 时 延 的 精度 也 可 以 得 到 


时 采用 该 


于 预测 大 气 对 深 空 站 链 路 中 Ka 波段 信号 的 影响 。 天 气 预报 数值 模型 用 于 预测 短期 大 气 状 
三 维 辐射 传输 模型 ， 用 于 计算 倾斜 路 径 衰减 和 天 线 噪 声 温 


| Ka 波段 的 接收 机 信 噪 比 。 相 比 传统 ITU-R， 该 方法 可 以 表征 热 不 均匀 和 大 气 


态 ; 微 物理 定向 传播 方案 主要 是 
RE, BET HUM 


定向 传播 方案 相 


结合 


的 方法 ， 用 


散射 。 利 用 区 域 天 气 预报 数值 模型 对 大 气 状 态 进行 短期 预测 ， 选 择 “ 贝 皮 科 伦 布 号 ”探测 器 
进行 测试 ， 在 接收 短期 和 丢失 数据 量 方面 ， 使 用 每 日 天 气 预 报 数值 模型 ， 而 不 是 每 月 的 气候 
统计 数据 ， 每 年 的 数据 量 回报 可 以 增加 2096 DLE”. 

2021 年 ， 考 虑 系统 性 测量 误差 及 模型 稳定 性 ， 北 京 航天 飞行 控制 中 心 对 中 性 大 气 天 项 
时 延 和 映射 函数 的 混合 模型 进行 了 改进 ， 并 在 “嫦娥 五 号 ”中 进行 了 验证 ， 有 具体 公式 如 下 : 
d 

TES RU i (8) 


Et, D 为 中 性 大 气 时 延 ，d、a 和 5b 为 模型 参数 ，9 是 观测 仰角 。 


H 


段 均衡 分 布 ， 提 高 了 低 仰角 弧 段 中 性 


另外 ， 国 际 地 球 自 
TERS) 规范 针对 无 线 电 和 光学 地 基 测 量 数据 的 


FP 性 大 气 时 延 混合 模型 解决 了 仰角 越 低 误 差 越 大 的 问题 ， 


FE 大 气 时 延 模型 精度 ; 佳木斯 、 
(HER RMS ( 含 实测 设备 误差 ) 最 大 差别 为 0.165 m、0.012 m^", 


使 得 建 模 误差 在 整个 观测 弧 
喀什 深 空 站 模型 值 与 实测 


转 服务 (International Earth Rotation and Reference Systems Service, 


i51] 


于 无 线 电 地 基 测 量 数据 ， 采 用 Davis 等 人 于 19 
值 修正 中 性 大 气 时 延 ， 然 后 采用 
VMFI 不 可 用 时 采用 GMF. Pe 
地 基 测 量 数据 ， 推 荐 采 


] Mendes 和 Pavlis 在 2 


型 与 射线 追踪 技术 获取 的 中 性 大 气 天 项 时 延 比 对 RM S 为 1 mm^. 


性 大 气 时 延 


E 荐 了 相应 的 修正 方法 


。 对 


85 年 提出 的 S 


VME 1 映射 函数 模型 将 天 顶 时 延 投 影 到 斜 路 径 方 向 ， 在 


aastamoninen 模型 作为 先 验 


生 大 气 湿 延迟 和 梯度 作为 未 知 
的 中 性 


004 年 提 昌 


参数 进行 估计 ”。 对 于 光学 
EX UR TRIS GERA 


型 ， 该 模 
然后 ， 采 用 Mendes 等 


3 期 周 伟 莉 ， 等 : 深 空 探测 中 地 基 测 量 数 据 传播 介质 误差 修正 方法 


人 在 2002 年 提出 的 映射 函数 模型 ， 将 天 项 时 延 投影 到 斜 路 径 方 向 ， 该 映射 函数 模型 基于 


Marini 连 分 式 表示 ， 采 用 180 个 全 球 探 空 站 1 a 的 序列 拟 合 模 型 系数 。 


2.3 小 结 


本 节 对 中 性 大 气 时 延 产 生机 理 和 修正 方法 ; 
的 优 缺 点 ， 表 1 给 出 了 几 种 常用 中 性 


比 对 情况 。 在 实时 任务 下 ， 可 以 采 
技术 进行 修正 ， 但 微波 辐射 计 的 缺点 是 下 月 
等 地 区 实测 资料 构建 ， 基 于 数值 预 
的 修正 精度 还 有 待 进一步 验证 。 


表 1 


及 模型 的 射线 追踪 技术 修正 
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行 了 介绍 。 为 了 更 直观 区 分 各 种 修正 方法 
EE 大 气 时 延 不 同 修正 方法 在 时 效 性 、 
j 拟 合 天 项 时 延 模型 + 映射 函数 、 微 波 辐 射 计 和 射线 追踪 
天 无 法 使 用 。 由 于 数值 预报 模型 大 部 分 


修正 精度 等 方面 的 


基于 欧美 


不 同 中 性 大 气 时 延 修正 方法 对 比 


FP 性 大 气 时 延 在 中 匡 


深 空 站 


修正 方法 


时 效 性 ”修正 精度 ”实测 资料 备注 


天 项 时 延 模型 十 映射 函数 


温 湿 压 气象 资料 十 映射 函数 
小 时 更 新 GNSS 十 映射 函数 


事后 GNSS 十 映射 函数 
微波 辐射 计 
射线 追踪 


实时 


事后 
实时 
实时 


n 
p 


低 
ESS 咯 { 
y 住 实时 rn 


不 需要 


需要 雨天 无 法 使 / 


不 需要 受 数值 预报 模型 精度 限制 


相 比 各 种 模型 ， 基 于 GNSS 技术 反 演 是 修正 中 性 大 气 天 顶 时 延 精 度 较 高 的 一 种 方法 。 
但 该 方法 缺点 是 还 不 能 用 于 实时 任务 。 目 前 GNSS 跟踪 站 可 以 提供 实时 数据 流 ，IGS 实时 


工作 组 可 以 提供 实时 轨道 钟 差 产 品 ， 基 于 实时 数据 流 获取 实时 中 性 
是 未 来 深 空 探测 中 地 基 测 量 数据 中 公 


目前 国际 上 GPT 系列 模型 修正 精度 相对 较 高 ， 可 以 应 用 于 VLBI 观测 实时 人 和 修正。 上 
的 映射 函数 模型 基于 全 球 数据 构建 ， 而 各 个 区 域 气象 条 件 差异 较 大 ， 为 了 进一步 提高 


FE 精 度 ， 可 以 考虑 开展 本 地 化 映射 函数 方面 的 深入 研究 。 


气 时 延 在 低 高 度 角 观 测 时 的 修 


3 电离 层 


本 文 所 讨论 电离 层 主要 指 


E Hh TR 


3.1 ”电离 层 时 延 产 生机 理 


电离 层 具有 非 均 匀 性 、 各 向 异己 


ERI E BU 


50 km ~ 1000 km 高 度 的 地 球 高 层 大 气 空 域 。 该 
存在 相当 多 的 自由 电子 和 离子 ， 使 无 线 电波 改变 传播 速度 ， 发 生 折 射 、 反 射 和 散射 等 。 


E 大 气 时 延 改 了 
大 气 天 项 时 延 改 正 的 一 个 研究 方向 。 另 外 ， 了 映射 函数 
的 精度 直接 影响 深 空 探测 中 地 面 站 低 高 度 角 观 测 时 地 基 测 量 数 据 中 性 大 气 时 延 修正 的 水 平 ， 


E 成 为 可 能 ， 


于 目前 


性 大 


区 域 


FE 特 征 。 非 均匀 性 表现 在 电离 层 


BT SEES 


空间 位 置 的 不 同 而 不 同 ,各 向 异性 
射 指数 依赖 于 波 的 频率 。 当 无 线 电波 穿 过 ! 
数 ， 而 是 光速 与 波 在 电离 层 中 传播 速度 的 比 


表现 为 折射 指数 随 波 的 传播 方向 不 同 而 不 同 ,色散 改 
包 离 层 时 ， 折 射 指数 不 再 等 同 于 真空 


值 。 相 比 真 空中 的 传播 ， 该 差异 称 为 | 


Exe IH UT 
FP 的 折射 指 
离 层 时 延 
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电磁 波 通过 电离 层 时 引起 的 传播 路 径 的 时 延 为 ; 


Ap- [nas [fas fal (9) 


其 中 ，ds 表示 信号 到 观测 台 站 的 真实 传播 路 程 ，dr 表示 信和 号 传播 的 几何 路 程 ， 第 二 个 因 式 
表示 几何 路 径 与 传播 路 径 间 的 弯曲 程度 大 小 。 

电离 层 中 带电 粒子 的 分 布 极 不 均匀 。 电 离 层 时 延 的 大 小 取决 于 航天 器 视线 、 观 测 频率 、 
一 天 中 的 时 间 、 季 节 和 自 当 前 太阳 周期 开始 以 来 的 时 间 。 在 S 波段 ， 天 项 电离 层 时 延 范围 
为 1 ~ 10 m， 在 高 度 角 约 0° 处 的 时 延 约 为 天 顶 时 延 的 3.5 fii. 

依赖 于 VLBI 技术 进行 探测 器 测控 时 ， 采 用 交 蔡 观测 探测 器 和 河 外 射电 源 的 方式 可 以 
消除 测 站 相关 的 误差 。 同 时 探测 器 与 河 外 射电 源 的 角 距 对 较 差 会 产生 影响 ， 即 角 距 越 大 ， 较 
差 误 差 越 大 。 电 离 层 时 延 对 较 差 产生 的 误差 可 通过 下 式 表示 : 


ATion RI v2 x am : Js . Osep , (10) 
C 


HP, cion 为 任意 方向 电离 层 时 延 修 正 后 的 误差 ，f, 是 映射 函数 的 导数 值 ，0,。, 为 探测 器 
相对 于 河 外 射电 源 的 角 距 。 
3.2 ”电离 层 时 延 修正 法 

电离 层 时 延误 差 与 观测 频率 有 关 ， 在 2 GHz 下 能 达到 几 十 纳 秒 。 电 离 层 时 延 与 观测 频 
率 的 平方 成 反比 ， 根 据 频率 不 同 常 分 为 双 频 观测 的 电离 层 时 延误 差 修 正方 法 和 单 频 观测 的 
电离 层 时 延误 差 修正 方法 。 
3.2.1 双 频 观测 电离 层 时 延误 差 修 正法 

假设 深 空 站 以 双 频 率 fa 和 fo 进行 探测 器 观测 ， 电 离 层 引起 的 无 线 电信 号 路 径 时 延 分 
别 为 : 


k 
Apion (fi) = B X CTE 
1 


k 
A Dion( f2) = j X Crp 


，Apion 为 电离 层 引 起 的 信号 时 延 ，ki = 40.28 m?-s 7, Cre 是 电磁 波 在 传播 路 径 上 的 


(11) 


kı x CTE 1 1 
on o ai S 


kı x Cre = Tiss 
c R- 


(t2 — tı) . (13) 
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然后 可 得 到 两 个 无 线 电 频 率 信 号 到 达 测 站 时 ， 由 于 电离 层 效 应 引起 的 时 延 为 ; 
Ation( 万 ) = f2 a (to — ti) 
1 ， 2 (14) 
f 
Ation(f2) = Pp? zx tı) 
TERS 规范 里 推荐 双 频 观测 时 采用 此 方法 修正 电离 层 时 延 的 影响 ”"。 对 于 深 空 探 测 地 基 


测量 数据 来 说 ， 可 以 通过 此 方法 消 
小 于 1mm， 可 以 忽略 。 美 
消除 了 地 球 电离 层 的 影响 
对 于 绝 大 部 分 深 空 探测 
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Ei 
IRZ MH 


VLBI 


国 “ 卡 西 尼 ”探测 器 第 


EA 


除 一 阶 电离 层 影 响 ， 残 余 的 | 


B A ES 


IRF 


灵 高 阶 项 对 群 时 延 的 影响 
一 个 采用 多 频 链 路 (Ka 和 X 波段 ) 观测 ， 


技术 定 轨 任 务 来 说 ， 所 用 8 


无 法 像 测 地 VLBI 技术 一 样 
不 适合 采用 双 频 法 实时 消除 1 
3.2.2 3 


= 


i TE 
AEE 


ET 


HZ 


对 无 线 电 信号 单 频 观测 来 说 ， 
与 频率 有 关 ， 电 离 层 时 延 修正 需要 区 

(1) 美国 探测 器 

1986 年 ，Lanyi 利用 GPS 卫 
型 ， 开 启 了 利用 GNSS 数据 反 演 
赛 ”轨道 确定 时 采 


e 


离 层 时 延 修正 法 
星 观 涡 
电离 层 模型 
j 深 空 站 ADOR 测量 


FILIM) E, A inc Ma 


电离 层 时 延误 差 可 以 通过 模型 
区 分 站 点 、 探 测 器 和 频率 。 


S 


上 数据， 采用 多 
的 新 方向 


时 数据， 电离 


层 改 正 模型 ; 行 修正 。 由 于 马 德 日 
测定 总 : 含量 5 


a 


号 探测 器 探测 金 和 
离 层 时 延 采 用 GPS 观测 数 


。2007 年 ，JPL 采用 GPS 
空 站 ADOR、VLBI 观测 值 和 VLBA BEAR 
时 采用 10 个 VLBA 天 线 跟 
居 反 演 的 垂直 总 电子 含 


通道 带宽 以 及 整个 频 
因 实 时 任务 无 法 保证 探测 器 观 测 信 标 全 部 开 


JPET BIE. 由于! 


PRAIA UM EE 
. 2001 Æ, ŽE 


启 ， 


IT 


带 跨度 
因此 并 


H 


层 时 延 


下 


EB SRDS EM 


观测 数据 反 演 


获取 的 上 


Uu 


Ma”. 


从 2006 年 到 2009 年 ， 美 


区 域 电离 层 改 了 
上 火星 探测 器 “ 奥 德 
EIN EAA GPS 观测 数据 反 演 的 电离 
有 E 深 空 站 没有 观测 数据 ， 采 用 差分 


Es 


技术 


层 改正 模型 标 校 深 


CE 


H» 


4/6 


Hoz 


DN 9 


全 球 分 布 图 进行 修了 
美国 深 空 探测 中 电离 层 时 延 修 


M 


其 他 频率 ， 计 算 每 个 测控 站 探测 器 
项 


~ 


E, VTEC 采样 间 


FERNS WE (AL 2295 MHz 频率 作为 基准 


JAR HE PY H 


ba 2h, X 


1E 


"Et 
视线 方 


E 务 需要 ,电离 层 校 准 H 


fH 
名 电离 层 校准 时 延 。 电 离 
illi 


式 表 示 ， 作 为 归 一 化 时 间 的 函数 ， 
须 提 供 探测 器 轨道 信息 。 
(2) 欧洲 探测 器 


欧洲 深 空 网 在 早期 主要 采 月 


HZ 


该 模型 基于 欧洲 定 轨 
全 球 总 电子 合 


一 2007 年 高 精度 4 


对 于 ADOR 测 
离 层 时 延 修正 法 
H Klobuchar 模型 


可 以 同时 提供 多 个 分 段 连续 多 : 


由 于 观测 不 包含 非常 ) 
量 (vertical total electron content, S EE 
方式 计算 测控 站 点 VTEC. 
针对 每 个 测控 站 点 、 航 天 器 和 频率 进行 设计 。 
对 延 可 以 很 容易 j 
层 校 准时 延 采 用 归 一 化 多 
项 式 。 在 使 


泛 分 离 的 频带 


IJ» 


已 


电离 层 校 
也 转换 到 


Jm, ] 


j 户 必 


电源 位 置 信 A oa 


量 数据 ， 还 必须 提供 河 外 射 ! 


4 修正 电离 层 时 延 “”， 
5096. 21 世纪 以 来 发 展 了 另 一 个 经 验 模型 NTCM-GL (Global Neustrelitz TEC Model): 
中 心 (Center for Orbit Determination in Europe, CODE) 提供 
含量 (total electron content, TEC) 数据 拟 合 而 成 ， 


该 模型 修 


正 精 度 为 


[13, sep 


tH 1998 
负 偏差 为 0.3 


TECU, RMS 偏差 为 7.5 TECU. 与 TOPEX/Poseidon 测 高 仪 结果 相 比 ，NTCM-GL 模 
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型 RMS 偏差 在 太阳 活动 低 年 和 高 年 分 别 为 5 TECU 和 11 TECU。 另 外 ， 基 于 合作 协议 ， 
欧洲 深 空 网 电离 层 时 延 修 正 也 由 美国 JPL 提供 技术 支持 。 

近 几 年 ， 除 了 中 性 大 气 时 延 修正 平台 ， 欧 洲 深 空 网 还 开发 了 一 个 独立 的 电离 层 时 延 修 
正平 台 。 该 平台 由 欧洲 导航 办 公 室 基于 GNSS 观测 数据 计算 电离 层 时 延 。 该 平台 尽量 依赖 
内 部 数据 和 软件 ， 简 单 、 可 靠 且 容易 操作 ， 使 用 NTCM-GL 模型 作为 背景 场 。 该 平台 需要 
的 输入 及 处 理 策略 包括 如 下 几 个 方面 : (1) 需要 欧洲 测控 站 的 GNSS 观测 数据 和 全 球 IGS 网 
观测 数据 ，(2) 需要 IGS 发 布 的 常规 GNSS 卫星 轨道 和 钟 差 产 品 、 常 规 斜 路 径 电 离 层 时 延 
(slant Total Electron Content, sTEC) 和 硬件 延迟 (Differential Code Bias, DCB); (3) 需要 
指定 48 h 间隔 的 欧洲 航天 器 /类 星体 轨迹 信息 ; (4) 欧洲 深 空 网 每 次 观测 不 会 超过 15 h, A 
此 每 天 的 电离 层 标 校 值 覆盖 48 h， 时 间 从 前 天 0 点 到 当天 0 点 ， 每 天 7 点 33 分 自动 运行 ; 
(5) 当 测 站 截止 高 度 角 > 5° 时 才 考 虑 进入 观测 ; (6) 同 中 性 大 气 标 校 一 样 ， 分 快速 模式 和 最 
终 模 式 ， 最 终 模 式 处 理 每 天 启动 时 间 与 快速 一 样 ， 区 别 在 于 延迟 4 d 启动 ; (7) 在 GNSS X 
测 数据 不 可 用 的 情况 下 ， 欧 洲 深 空 网 使 用 NTCM-GL 计算 电离 层 时 延 。 

GNSS 接收 机 和 卫星 DCB 从 导航 办 公 室 每 天 运行 的 全 球 电离 层 模型 估计 中 获得 ，DCB 
不 能 在 绝对 意义 上 进行 估计 ， 通 过 指定 卫星 DCB 总 和 为 零 的 基准 实现 。 为 了 避免 DCB fii 
计 值 出 现 显著 的 逐日 变化 ， 例 如 由 于 某 一 天 的 卫星 丢失 而 改变 平均 值 ， 将 过 去 7 d 的 每 日 
DCB 估计 值 合并 。 


由 于 探测 器 观测 频率 不 同 ， 通 过 以 下 公式 在 不 同 频 率 间 进行 电离 层 时 延 转换 。 
+40.3 
Jion = EH : SCTE ; (15) 


其 中 ，f 是 频率 ，gios 是 f 频率 下 的 电离 层 时 延 ， 在 计算 时 会 设置 f = 2295 x 10° Hz. 

另外 ， 欧 洲 “Smart-1” 探 测 阶段 、 美 国 火 星 探测 器 “Rover B” KME RAE” 
探测 器 在 土星 上 的 位 置 等 采用 EVN 进行 跟踪 。 在 VLBI 技术 数据 处 理 中 ， 采 用 IGS 提供 的 
采样 率 2 h ffi VTEC 格 网 文件 和 单 层 投影 函数 修正 电离 层 时 延误 差 ， 该 产品 Curg 精度 为 
2~9 TECU"。 具 体 计算 公式 如 下 : 


CvTE 
一 1 
Cre cos(z’) ' (16) 


其 中 ，Czs 为 观测 方向 的 倾斜 总 电子 含量 ，Cvra 为 垂直 总 电子 含量 ，z' 是 电离 层 单 层 模型 
到 目标 探测 器 的 天 项 距 。 


R 
z’ = arcsin (z rà z) -sin(z) , (17) 


RF, Ree, AH 为 电离 层 的 高 度 ，z 是 目标 探测 器 的 视 真 空 天 项 距 。 
每 个 测 站 上 电离 层 时 延 可 以 通过 下 式 计 算得 到 : 
5.308 018 Cre 


Tion = Ax P , (18) 


测 消除 一 阶 电 离 层 时 延 影响 。 当 探测 器 频率 发 生变 化 不 能 采 
延 时 采用 CODE 提供 的 全 球 ， 
据 反 演 构建 ， 采 


数 


3 期 


AER, S: 深 空 探 测 中 地 基 测 量 数据 传播 介 


(3) 日 本 探测 器 上 
2001 年 以 来 ， 


e 


H 


(4) 我 国 探测 


北京 航天 飞行 控 各 


离 层 时 延 修 正法 
本 采用 多 频 VLBI 技术 对 SELENE RI 


质 误差 修正 方法 


1 器 进行 轨道 确定 ， 通 i 
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过 双 频 观 


ar E 


涉 测量 需求 ， 结 合 深 空 


测 
T 


PR 
H 
整 ; 


踪 弧 段 内 目标 高 度 角 等 信息 ， 经 
中 继 星 任务 期 间 的 实测 数 所 
弧 段 ， 该 方法 优势 更 为 明显 
由 于 IGS ERER PI ZS TE uS] > 
态 网 ”加 密 的 方式 构建 
驻 域 稀 玻 的 问题 ， 提 高 了 我 
” FR AGE |B 
电离 层 时 延 后 ，CVN 站 6 
实测 数据 分 析 ， 与 CODE 的 全 球 电离 层 模型 产品 i 
定 轨 残 差 在 SI 波段 由 6.07 ns 降 
段 由 0.70 ns 降低 到 0.60 ns”. 
另外 ，IERS 规范 也 
对 于 事后 用 户 可 以 采用 IGS 提供 的 全 
CODE 提供 的 全 


中 
24 
修 


验 


响 。 
于 实时 用 户 可 以 采用 


zi 


PI 
比 
效 


控 模 式 实时 
跟踪 弧 段 前 


数 型 投影 函数 ; 
特性 ， 以 日 积 和 


(4) 依据 模型 系数 及 跟 


BEJE Cre 模型 修正 。 
115 阶 次 球 谐 系数 形式 表示 ， 时 间 分 辨 率 2 h， 通 
1 min 采样 间隔 电离 层 修正 数据 。 
离 层 时 延 修 正法 
1 中心 针 对 “嫦娥 四 号 ”中 继 星 任务 S 频段 信 标 信号 的 高 精度 实时 干 


该 


j 双 频 组 合 消 除 一 阶 电离 
电离 层 模型 基于 全 球 GNSS 网 观测 
过 平滑 内 插 方式 得 到 


层 时 ， 


空 测控 干涉 测 


量 系统 采 
干涉 测量 的 电离 层 时 延 修 正方 法 。 
3 d 左右 内 的 深 空 站 Cre 实测 数据 ; 
fti] TEC 为 视线 向 的 测量 值 ， 首 先 依 和 
(3) WES 


E E R 
EYJA Cyrs 进行 质量 评估 
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该 方法 
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在 5096 AA”. 
3.23 ”电离 层 时 延误 差 修 正经 验 模型 


地 基 测 量 数据 中 1 
考 电 离 层 (international reference ionosphere, IRI) 模型 、 
验 模 型 的 优点 是 


M 经 验 模 型 


51,63] 


ELS ENTE Et 


. 2 


CVN 站 上 电离 层 时 延 修正 精度 。 医 


CT S H | CODE 和 SHAO 提供 


P, S1 波段 采 


电离 层 精 修正 模型 SHAO-GIM， 该 方法 解决 了 IGS 跟踪 网 在 


的 稀 玻 标 校 工作 模式 ， 研 究 验证 了 一 种 面向 
如 下 几 个 方面 
由 于 测 站 配置 的 GPS 双 频 接收 机 提 
影 函 数 转化 为 天 项 向 Cvrn 值 ， 选 取 三 角 
， 剔 除 异 党 时 间 段 观测 量 ， 基 于 周 


(1) 基 


UOS ELSE. Curs 为 变量 拟 侣 多 项 式 模型 ， 模 型 阶 数 可 视 拟 合 残 差 动态 调 
院外 段 的 日 积 秒 ， 计 算 测 让 Cvrs ll. (5) 结合 数据 处 理 系统 中 四 
影 函 数 转化 计算 视线 向 的 Ca, 预测 值 。 通 过 “嫦娥 四 号 
分 析 ， 该 方法 可 以 将 实时 测量 精度 提升 ~ nes DRE 
p, [61] 
跟踪 站 间距 大 。 中 国 科学 院 上 海天 文 台 CVN 采用 中 国 “ 陆 


4 是 2014 年 10 H 


的 GIM 模型 


6 条 基线 方向 的 定 轨 残 差 。 采 
氏 到 3.33 ns; 


荐 采用 Cvrs 模型 和 映射 函数 组 合 的 方式 修 


Ek VTEC 预报 模型 ， 全 球 Cure 


] €€ 嫦娥 五 号 [= ?9 
行 对 比 ， 
在 S2 波段 由 


HARE 
基于 SHAO-GIM 的 VLBI 
9.10 ns 降低 到 7.07 ns; 


ERR 


XO 3X 
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球 格 网 Cyrn 模型 ， 


正 电 离 
TERRE 10% ~ 20%. X 
E 预报 模型 平均 修正 误 


层 时 延 的 影 
对 


如 ，IRI 模型 
Bo 5095, 


融合 了 1 
另外 ，3.1 节 提 到 的 


电离 


Klobuchar 模型 、 


一 些 经 验 模型 予以 修正 。IERS 规范 推荐 采用 国际 参 
NeQnuick 模型 或 


E3 


层 探 测 仪 、 非 相干 散射 雷达 、 卫 星 资 料 和 探 空 火箭 资料 等 ， 
欧洲 深 空 探测 中 采用 的 背景 场 模 型 NTCM-GL 


可 以 应 用 于 深 空 探 测 实时 任务 ， 缺 点 是 修正 精度 低 。 


改正 
司 样 是 
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ik: 红色 表示 CODE; 蓝 色 表示 SHAO. 
4 2014 年 10 月 24 日 “ 媒 娥 五 号 ”再 入 返回 飞行 试验 中 ，S1 波段 采用 CODE 和 SHAO 模型 修正 电 
离 层 时 延 后 CVN 站 6 条 基线 方向 的 定 轨 残 差 


一 种 经 验 模型 。 
3.3 ”小 结 
本 节 介 绍 了 电离 层 时 延 产生 机 理 和 目前 常用 的 几 种 修正 方法 。 为 了 更 直观 地 区 分 常用 
电离 层 时 延 修 正方 法 之 间 的 差异 ， 表 2 在 时 效 性 、 修 正 精 度 和 观测 频率 要 求 等 方面 进行 了 
对 比 。 对 于 实时 单 频 观 测 ， 建 议 采 用 经 验 模型 或 全 球 VTEC 预报 模型 + 映射 函数 修正 电离 
层 时 延 ， 对 于 事后 任务 ， 推 荐 采用 全 球 VTEC 事后 模型 + 映射 函数 组 合 的 方式 修正 电离 层 
时 延 ， 对 于 有 双 频 观测 的 次 空 探测 任务 ， 可 以 采用 双 频 组 合法 修正 电离 层 时 延 的 影 


a 


o 


表 2 常用 电离 层 时 延 修正 方法 对 比 


修正 方法 时 效 性 ”修正 精度 。 观测 频率 要 求 
双 频 组 合法 实时 ”消除 一 阶 项 WUT DE 
全 球 VTEC 事后 模型 十 映射 函数 ”事后 高 单 频 及 以 上 
ESR VTEC 预报 模型 十 映射 函数 ”实时 低 单 频 及 以 上 
Klobuchar 模型 实时 低 单 频 及 以 上 
NTCM-GL 模型 实时 低 单 频 及 以 上 
NeQuick 模型 实时 低 单 频 及 以 上 
IRI 模型 实时 低 单 频 及 以 上 
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H 


近 几 年 来 ，IGS 开始 提供 全 球 实时 电离 层 产 品 ， 采 用 15 阶 次 球 谐 函数 形式 ， 时 间 分 辨 
率 为 20 min。 与 最 终 全 球 电离 层 模型 产品 比较 ， 精 度 约 为 3 TECU ”， 是 目前 电离 层 时 延 
修正 研究 的 一 个 热点 。 随 着 电离 层 研究 的 不 断 深入 ， 小 时 更 新 及 实时 全 球 电离 层 模型 产品 得 
到 发 展 并 有 望 应 用 到 深 空 探测 任务 。 该 类 产品 可 以 提高 电离 层 时 延 修正 的 精度 及 时 效 性 。 


4 太阳 等 离子 体 


探测 器 在 行星 际 空间 飞行 时 ， 会 受到 太阳 风 的 强烈 影响 。 太 阳 风 是 高 密度 强 汝 流 电 离 气 
体 (等 离子 体 ) 从 太阳 中 喷射 出 来 的 结果 ， 这 些 电离 的 离子 以 400 km/s 的 速度 从 太阳 高 速 流 
出 后 扩散 到 广阔 的 行星 际 空间 ， 使 行星 际 空间 成 为 等 离子 体 区 ""。 当 电磁 波 穿 过 这 些 区 域 
时 ， 通 信 信 号 产生 时 间 延 迟 ， 最 大 可 达 几 十 米 。 

太阳 等 离子 体 时 延误 差 通常 采用 经 验 模 型 估计 。1981 年 ，Muhleman 和 Anderson 提出 
了 M&A 电子 密度 模型 "" ， 该 模型 可 以 较 好 地 反映 等 离子 区 物理 背景 ， 在 工程 任务 中 应 用 
较 多 。 美 国航 空 航天 局 提供 的 日 冕 距离 修正 公式 使 用 了 M&A 模型 。 该 模型 中 等 离子 体 密 度 
可 近似 表述 为 : 


D p? 
Ne(r, p) = rT exp ( Z) | 2 , (19) 


" 
其 中 , r 为 日 心 到 探测 器 与 望远镜 连 线 的 径 向 距离 (单位 为 太阳 半径 ); v 为 该 点 的 日 面 纬 
度 ， 即 某 点 与 日 心 的 连 线 和 太阳 赤道 面 的 夹 角 ;ws 为 日 心 距离 相同 处 等 离子 体 密度 最 大 值 
所 对 应 的 日 面 纬 度 : D B 为 模型 系数 ， 由 观测 数据 拟 合 得 到 。 Muhleman 和 Anderson 
利用 火星 上 的 Viking 飞船 的 测量 数据 得 出 了 1976 年 太阳 极 小 期 的 电子 密度 ， 他 们 取 
D — 1.32 x 106, B = 2.3 x 10°. 

1995 4E, Bird 4& A" FEF Ulysses 飞 船 的 双 频 信号 进行 反 演 并 建立 了 Bird 模型 。2010 
年 ， 加 州 理 工学 院 研究 人 员 提 出 Caltech 模型 "" ， 该 模型 适用 于 低 日 面 纬度 地 区 ， 只 与 目标 
到 日 心 的 距离 有 关 ， 运 算 逻 辑 简单 。 该 模型 公式 如 下 : 


—6 一 2.3 
No) = 221 x 10! (£) +155xao (总) 20 
(r) [3 T, (20) 


其 中 ，Ne(r) 为 等 离子 体 密度 ，a 为 日 心 到 探测 器 与 望远镜 连 线 的 径 向 距离 (单位 m)，Ro 
为 太阳 半径 (Ro = 6.96 x 10° m)。 图 5 是 采用 Caltech 模型 计算 太阳 等 离子 体 时 的 几何 关系 


刘 庆 会 等 人 ”对 采用 VLBI 技术 进行 “天 问 一 号 ”测定 轨 时 等 离子 体 时 延 的 影响 进行 
了 仿真 分 析 ， 利 用 “天 问 一 号 ”探测 器 发 射 到 着 陆 过 程 中 太阳 -地 球 -探测 器 夹 角 的 变化 进行 
计算 ， 发 现 太阳 等 离子 体 时 延 相对 较 小 。 北 京 航天 飞行 控制 中 心 采 用 Caltech 模型 和 M&A 
模型 分 别 对 “天 问 一 号 ”太阳 等 离子 体 时 延误 差 进行 了 分 析 ， 结 果 表 明 两 种 模型 精度 相当 。 
图 6 是 “天 问 一 号 ”探测 时 佳木斯 、 喀 什 和 阿根廷 深 空 站 等 离子 体 延迟 对 UXB 测 距 的 双 程 
影响 误差 ， 最 大 约 18.6 m。 在 采用 测 距 测速 数据 进行 定 轨 时 需 考虑 等 离子 体 影响 ” 。 
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图 5 ”计算 太阳 等 离子 体 的 几何 关系 


óp./m 


YE: 蓝 色 表 示 佳 木 斯 深 空 站 ， 红 色 表 示 喀 什 深 空 站 ， 绿 色 表 示 阿 根 廷 深 空 站 。 


Elo ”等 离子 体 延 迟 对 UXB 测 距 的 双 程 影响 误差 
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or 


Ej, Se: 深 空 探测 中 地 基 测 量 数据 传播 介质 误差 修正 方法 


深 空 探测 中 探测 器 轨道 确定 时 ， 传 播 介质 误差 是 一 个 很 重要 的 误差 源 ， 主 要 包括 中 性 大 
气 时 延误 差 和 电离 层 时 延误 差 等 。 为 了 有 效 削弱 该 误差 源 的 影响 ， 根 据 时 效 性 差异 ， 需 要 选 
性 大 气 和 电离 层 时 延误 
则 控 站 ， 俄 罗斯 近 三 十 


Z 


择 合适 的 中 性 大 气 和 电离 层 时 延 改正 方法 。 为 此 ， 我 们 详细 分 析 了 中 
差 产生 机 理 及 各 测控 站 修正 方法 。 主 要 包括 美国 、 欧 盟 、 日 本 和 我 国 


uu 
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年 探测 任务 相对 较 少 ， 且 公开 资料 甚 少 。 对 于 中 性 大 气 时 延误 差 来 说 ， 通 常 可 以 采用 天 项 时 
延 十 映射 函数 模型 和 基于 WVR 修正 法 等 进行 改正 。 若 探测 器 有 双 频 信和 号 观测 ， 可 以 采用 双 


频 组 合 消除 一 般 电 离 层 影响 ， 高 阶 项 影响 很 小 ， 可 以 忽略 ; 然而 绝 大 部 分 探测 器 采用 单 频 信 
号 观测 ， 此 时 可 以 采用 全 球 或 区 域 电离 层 模型 进行 修正 。 此 外 ， 对 于 行星 际 的 探测 ， 会 受到 


等 离子 体 影 响 ， 我 们 简要 介绍 了 等 离子 体 误差 的 常用 修正 方法 。 
随 着 深 空 探测 的 不 断 深 入 ， 对 探测 器 定 轨 精度 的 要 求 也 越 来 越 高 。 由 于 中 | 


Py 


EUN GERE 


响 ， 低 高 度 角 定 轨 时 延 残 差 仍 然 较 大 ， 是 制约 探测 器 定 轨 精度 提升 的 一 大 因素 。 目 前 采用 的 
都 是 全 球 通 用 的 映射 函数 模型 ， 但 各 个 测 站 气象 条 件 差异 较 大 。 所 以 ， 针 对 深 空 站 的 本 地 化 


映射 函数 模型 的 构建 是 未 来 的 研究 重点 ， 也 是 国际 上 的 研究 热点 及 难点 。 
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An Review of the Propagation Delay Calibration of Ground 


Station Observations in the Deep Space Exploration 


ZHOU Wei-li-?3, SONG Shu-li?, HU Xiong!, JIANG Jun?’ 


(1. National Space Science Center, Chinese Academy of Sciences, Beijing 100190, China; 2. Shanghai 
Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China; 3. University of 
Chinese Academy of Sciences, Beijing 100049, China) 


Abstract: In deep space exploration, the spacecraft trajectory is mainly measured by the 
ground tracking stations. The radio signals transmission between spacecraft and stations will 
be affected by the earth's propagation medium error, which contains neutral atmospheric 
delay and ionospheric delay, et al. According to the application requirements of propagation 
delay calibration in deep space exploration, we investigate the ground station observations' 
propagation delay at home and abroad, and focus on the basic principles and calibrated 
methods of VLBI observations! neutral atmosphere and ionosphere. The development trend 
of the key technologies in the future are analyzed, to improve the accuracy and timeliness of 
the propagation delay calibration. This paper can provide a reference for propagation delay 


calibration for trajectory determination of the spacecraft. 


Key words: deep space exploration; neutral atmospheric delay; ionospheric delay 


